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Abstract — Super-spreading events for infectious dis- 
eases occur when some infected individuals infect 
more than the average number of secondary cases. 
Several super-spreading individuals have been identi- 
fied for the 2003 outbreak of severe acute respiratory 
syndrome (SARS). We develop a model for super- 
spreading events of infectious diseases, which is based 
on the outbreak of SARS. Using this model we de- 
scribe two methods for estimating the parameters of 
the model, which we demonstrate with the small-scale 
SARS outbreak at the Amoy Gardens, Hong Kong, 
and the large-scale outbreak in the entire Hong Kong 
Special Administrative Region. One method is based 
on parameters calculated for the classical suscepti- 
ble - infected - removed (SIR) disease model. The 
second is based on parameter estimates found in the 
literature. Using the parameters calculated for the 
SIR model, our model predicts an outcome similar to 
that for the SIR model. On the other hand, using 
parameter estimates from SARS literature our model 
predicts a much more serious epidemic. 

Keywords: modeling, infectious diseases, super- 
spreading events, severe acute respiratory syndrome 

1 Introduction 

Severe acute respiratory syndrome (SARS) is a highly 
contagious respiratory disease which is caused by the 
SARS Coronavirus. It is a serious form of pneumonia, re- 
sulting in acute respiratory distress and sometimes death. 
SARS emerged in China late 2002 and quickly spread 
to 32 countries causing more than 774 deaths and 8098 
infections worldwide. One of the intriguing characteris- 
tics of the SARS epidemic was the occurrence of super- 
spreading events. Super-spreading events for a specific 
infectious disease occur when certain infected individu- 
als, called super-spreaders, infect more than the average 
number of secondary cases. 

According to the U. S. Centers for Disease Control and 
Prevention, a person is a super-spreader if they cause 
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more than 10 secondary infections. Such super-spreading 
individuals have been identified in the SARS outbreak 
and they are thought to have caused most of the sec- 
ondary infections. For example, in Singapore, about 
80% of infections have been attributed to only 5 super- 
spreading individuals ([7]). One extreme super-spreading 
individual in Hong Kong caused more than 100 secondary 
infections ([13). To contrast, in Singapore, most individ- 
uals caused secondary infections ([10]). 

One differential equation model for super-spreading indi- 
viduals was proposed before the recent SARS outbreak, 
by Kemper ([BJ). He presents a modified susceptible - 
infected - removed (SIR) disease model to capture the 
effect of the super-spreading individuals in which the 
infected individuals are split into two different infected 
classes with different transmission rates. The super- 
spreaders have a higher transmission rate, meaning more 
of their contacts with susceptible individuals result in a 
new infection than the regular infected individuals. This 
was the first model designed specifically to address the 
effect of super-spreading individuals on the course of in- 
fectious disease epidemics. 

Due to the extreme influence of super-spreading individ- 
uals during the SARS outbreak, many mathematical epi- 
demiologists developed models for the spread of SARS, 
which included the occurrence of the super-spreading 
eventqj. It is interesting to note how each incorpo- 
rated the super-spreading events into their model. First 
are models where the super-spreading individuals have a 
higher transmission rate than the regular infected indi- 
viduals, as in Kemper ([BJ). For example, Masuda, et 
al. ([ID]) developed a contact network model which had 
different transmission rates for regular infected individu- 
als and for super-spreading individuals. Second are mod- 
els in which parameters are taken from probability dis- 
tributions, making the super-spreading individuals nat- 
urally appear as the right-hand tail of the distributions. 
Lloyd-Smith, et al. ([9]) developed a stochastic compart- 
ment model and Fang, et al. ([3]) developed a spacial lat- 
tice combined with a deterministic compartment model 
wherein the individual reproduction numbers are drawn 
from a continuous probability distribution. 



x Many models were developed to try and mathematically cap- 
ture the SARS outbreak. Interested readers are directed to Bauch, 
et al. (ff]) for an overview of SARS models. 



The goal of this paper is to capture the effect of the super- 
spreading individuals using a modification of the classi- 
cal SIR disease model. The modification is inspired by 
Li, et al. (|S]) who determined that for SARS "the daily 
infection rate did not correlate with the daily total num- 
ber of symptomatic cases but with the daily number of 
symptomatic cases who were not admitted to a hospital 
within 4 days of onset of symptoms." This means that the 
number of infected individuals is closely associated with 
the number of individuals who remain outside of isolation 
longer than most other infected individuals. These indi- 
viduals remaining outside of isolation longer than normal 
are related to the disease severity. They have more con- 
tacts with susceptible individuals and thus more chances 
to spread the disease, in other words, they are the super- 
spreaders. In our model, we split the infected individuals 
into two classes, with two different "removal" rates; these 
two rates determine how long an infected individual re- 
mains outside of isolation. 

In this paper, we modify the classical SIR disease model 
to capture the effect of a super-spreading event, using the 
idea that super-spreading individuals stay out of isolation 
longer than individuals who are not super-spreading. A 
description of the SI PR model, the model equations, and 
some basic model properties are given in Section [5] We 
demonstrate the model using a small-scale and a large- 
scale outbreak of SARS in Section [3] Our discussion fol- 
lows in Section |4j 



1.1 The SIR model 

We conclude the introduction with a brief review of the 
classical SIR disease model. A more detailed review is 
given in Murray ([H]), and in Ching, et al. ([2]) and the 
references given there. 

The classical SIR disease model splits a fixed-size popu- 
lation, N, into three distinct classes: the susceptible indi- 
viduals, S, those do not have the disease and can become 
infected; the infected individuals, /, those who have the 
disease and can infect susceptible individuals; and the re- 
moved individuals, R, those who have recovered, die, or 
moved into isolation. Individuals in the removed class 
gain permanent immunity and remain in the R class for- 
ever. Schematically, the SIR model is 




Figure 1: SI PR model schematic, where S are the sus- 
ceptible individuals, I are the regular infected individu- 
als, P are the super-spreaders, R are the removed indi- 
viduals, /3 is the transmission rate, b is the probability 
that a new infection will be a regular infected person, 
1 — b is the probability that a new infection will be a 
super-spreading individual, v\ is the removal rate for a 
regular infected individual, and vi is the removal rate for 
a super-spreading individual. 



equations describing the SIR model is 
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where /3 is the transmission rate and v is the removao 
rate. 

2 The SIPR model 

We describe the SIPR model, a modification of the clas- 
sical SIR model, which captures the effect of super- 
spreading individuals. Schematically, the model is given 
in Figure [TJ While the basic SIR model has one class of 
infected individuals, the SIPR model has a second class 
of infected individuals, the super-spreaders, denoted by 
the variable P. In the SIPR model we divide the popu- 
lation, of fixed size N, into four groups, namely the sus- 
ceptible individuals, S, the regular infected individuals, 
I, the super-spreaders, P, and the removed individuals, 
R. 



Individuals in the population are assumed to homoge- 
neously mix. Contacts between the susceptible and in- 
fected individuals result in a new infected individual at a 
rate proportional to the to the number of susceptible and 
infected individuals. Infected individuals are removed to 
class R at a rate proportional to the number of infected 
individuals. The system of nonlinear ordinary differential 



A susceptible individual can become infected through 
contact with either a regular infected individual or a 
super-spreading individual. Then with probability b the 



2 The R class is often called the recovered class containing those 
who have recovered from the disease. In this situation, all indi- 
viduals recover from the disease; no individuals die or move into 
isolation. The parameter v is correspondingly called the recovery 
rate. 



new infected individuals become regular infected individ- 
uals (move to class /) and with probability 1 — b the new 
infected individuals become super-spreading individuals 
(move to class P). We assume most new infections are 
regular infected individuals, that is, b > 1 — 6. We as- 
sume that the two infected classes, I and P, have the 
same transmission rate, /3. 

To capture the effect of the super-spreading individuals, 
we use the idea that individuals who are super-spreaders 
stay out of isolation longer than the regular infected 
individuals. The regular infected individuals stay out 
of isolation (the R class) for l/v\ days and the super- 
spreading individuals stay out of isolation for 1/j/2 days, 
with l/v\ < I/2/2. Therefore, the model has two distinct 
removal rates namely, v\ and v 2 , corresponding to regu- 
larly infected individuals, /, and super-spreading events, 
P, respectively. When a person is removed to the R class 
there is no possibility of becoming susceptible again, but 
rather they recover and gain permanent immunity, or die; 
in either case, they remain in the removed compartment 
forever. 

Based on the previous descriptions and assumptions we 
formulate a system of four ordinary differential equations 
for the SI PR model. 



dS 

~dt 
dl 
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dP 
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-P{I + P)S 

bj3(I + P)S - vil (1) 

(1- b)0(I + P)S-u 2 P 

V\I + v 2 P 

subject to the following initial conditions 
5(0) = N - I - P = S Q , 

1(0) = I , P(0)=P o , and P(0) = 0. 

The SIPR model given by Equations ([1} has the follow- 
ing properties. 

• The SIPR model has a unique global solution. 

• The components of the solution, S(t), I(t),P(t), 
and R(t) , of the SIPR model are non- negative and 
bounded by N for all time, t > 0. 

• The SIPR model has equilibrium points (N, 0, 0, 0), 
{S*,0,0,R*), for any 0< S* < S , and (0,0,0,7V). 

• The individual reproduction number is 
bf3S , (l-b)0S o 



Ro = 



V\ 



v-i 



The properties of the SIPR are, in general, the natu- 
ral modification of the corresponding properties of the 



classical SIR model. It is interesting to note that the 
Ro for the SIPR model is the Rq for the corresponding 
SIR model for each of the two infected classes, I and P, 
multiplied by the probability of a new infected individual 
becoming an I or a P, respectively. 

The SIPR model can be used to analyze any infectious 
disease where super-spreading events have been identi- 
fied. As a particular example, super-spreading events 
have been identified in one outbreak of measles as de- 
scribed in Paunio, et al. ([T2]). 

3 The SARS Epidemic 

We use the SIPR model to study the spread of SARS on 
a small scale, in the Amoy Gardens apartment complex 
in Hong Kong, and on a large scale, in the entire Hong 
Kong Special Administrative Region. In both cases, we 
fit the model to the data using two parameter estimation 
methods. We begin with a description of the estimation 
procedures. 

The general parameter estimation procedure was as fol- 
lows. The SIPR system of equations ([1]) was solved re- 
peatedly with parameter sets taken from the allowable 
range of possible parameters, and the least squares error 
between the cumulative number of cases of the solution 
and the actual data was computed. The least squares 
error was minimized. All computations were done with 
MATLAB using the function fminsearchbnd. 

In parameter estimation Method 1, we used the param- 
eters estimated with a fit of the classical SIR model to 
the data. For the data on the entire Hong Kong Special 
Administrative Region, we were unable to find SIR pa- 
rameters. Therefore, we began with fitting the classical 
SIR model to the data. A summary of the SIR model 
parameters is given in Table [T] We use the transmission 
rate of the SIR model as the transmission rate of the 
SIPR model. The transmission rate, the removal rate, 
and the initial number of susceptible individuals are used 
to determine the basic reproduction number P ■ We esti- 
mate parameters b, v\ , and v 2 so that together they have 
Ro as calculated, and they satisfy b > 1 — b and v\ > v 2 . 

In parameter estimation Method 2, bounds for the pa- 
rameters in the SIPR model were determined from the 
literature. Then, using these bounds, we estimate the 
parameters that give the best fit to the data. Details and 
references for these bounds are given in the next para- 
graph and in Subsections 13.11 and 13.21 

For estimation Method 2, we make some assumptions re- 
garding the parameters that apply to both the small- and 
large-scale outbreaks. In situations with super-spreading 
individuals, it is assumed that most infected individuals 
are not super-spreaders, which is the case for the spread 
of SARS. In fact, Masuda, et al. ([ID]) state that for the 



Parameter 



Amoy Gardens Hong Kong Hong Kong 

Feb 21 - Mar 21 Mar 29 - Jun 12 



v 

Ro 



1.4850/N 
0.9750 
1.5154 



0.2586/N 
0.0821 
3.1511 



0.1351/N 
0.1923 
0.7025 



Table 1: Estimated parameters for the SIR model and the SARS outbreak in the Amoy Gardens, Hong Kong (from 
([2])), and in the Hong Kong Special Administrative Region for both February 21 - March 21, 2003 and March 29 - 
June 12, 2003. 



SARS outbreak in Singapore 80% of infected individuals 
infected no one else. This leads us to set a lower bound on 
the probability of becoming a regular infected individual 
of 0.8, and we have 0.8 < b < 1. For SARS, after on- 
set of symptoms regular infected individuals stayed out 
of isolation between 3 and 5 days ([3]); we assume the 
average and set V\ = 1/4. We assume that on average 
super-spreading individuals moved into isolation 10 days 
after they became infectious, that is, v% = 1/10. Finally, 
many researcher have shown that the individual repro- 
duction number, Ro, for the SARS outbreak is between 
1.5 and 4 ([IJ). We assume this region for Rq, which gives 
a constraint on the parameters /3, b, i>i, and v^. 



3.1 The SARS epidemic - small scale 

We use the SIPR model to study the spread of SARS at 
the high rise apartment building Amoy Gardens, Block E, 
in Hong Kong. We assume that the cumulative number of 
confirmed cases of SARS in the Amoy Gardens, summa- 
rized in Table [2] are all from Block E, and that Block E 
contains a total of 792 individuals, as in Ching, et al. ([2]). 

For parameter estimation Method 1, we use the param- 
eters determined by Ching, et al. ([2]). The fixed and 
estimated parameters of the SIPR model are given in Ta- 
ble |31 The resulting number of confirmed cases of SARS 
are presented in Table [2] 

For parameter estimation Method 2, we survey litera- 
ture to determine appropriate bounds on all of the SIPR 
model equations. We assume that the total number of in- 
fected individuals at t = is 4, as in Ching, et al. ([2]). 
We conservatively assume < j3 < 0.01. The bounds 
and constraints for each parameter as well as the esti- 
mated parameter values are given in Table [3] The result- 
ing number of confirmed cases of SARS are presented in 
Tabled 



3.2 SARS epidemic - large scale 

We use the SIPR model to study the spread of SARS in 
the entire Hong Kong Special Administrative Region. We 
use the cumulative number of confirmed cases of SARS 
in Hong Kong as reported by the World Health Organi- 



zation^ (WHO). The first case of SARS in Hong Kong 
appeared on February 21, 2003. The WHO began daily 
reporting of SARS cases on March 17, 2003, and June 12, 
2003 shows the last additional confirmed case to a total 
of 1755 confirmed cases. In the parameter estimations, 
we use all of this data, and we summarize the data in 
Table H 

To begin, we identify three distinct time periods of the 
SARS epidemic in Hong Kong. The first time period is 
February 21 - March 21, 2003, the second is March 22 
- March 28, 2003, and the third is March 29 - June 12, 
2003. The first reported case of SARS in Hong Kong oc- 
curred on February 21, 2003, when a SARS infected doc- 
tor arrived from China. From this date until March 21, 
2003, very little was known regarding SARS and few pre- 
cautions were taken by individuals to prevent the spread 
of infection. On March 21, Hong Kong officials decided 
to begin daily reports to the public regarding the threat 
of SARS infection. As a result of these warnings people 
began to protect themselves from infection and the dy- 
namics of the disease transmission changed. Due to the 
continued spread of SARS and the increasing awareness of 
the seriousness of infection, on March 29, all Hong Kong 
schools were suspended, again changing the disease dy- 
namics. Classes resumed for all students by May 19 and 
the last reported case of SARS in Hong Kong occurred 
on June 12, 200fl 

We examine the first and third periods only, using the 
SIPR model. The first time period corresponds to quick 
spreading of the disease and we show this region has an 
Rq value above 1. The third period corresponds to the 
end of the disease and correspondingly has an Rq value 
below 1. 

In 2003, the Hong Kong area had a total population of 
6.803 million people. The number of susceptible individ- 
uals in Hong Kong during the SARS epidemic is not the 
entire population; the number of susceptible individuals 
must be approximated. Katriel and Stone ([5]) give a 



3 The WHO numbers were accessed on September 29, 2010, on- 
line at http://www.who.int/csr/sars/country/en/indcx.html 

4 A chronology of the SARS epidemic in Hong Kong can be found 
in Appendix III of the final report by the SARS Expert Committee 
of the Hong Kong Special Administrative Region, accessed October 
15, 2010, on-line at http://www.sars-expertcom.gov.hk/ 



Day 



March 26 March 27 March 28 March 29 March 30 



Confirmed 
Method 1 
Method 2 



7 
12.2 
9.1 



22 
26.2 
19.1 



56 
48.2 
37.4 



78 


114 


80.3 


123.8 


69.7 


122.3 



Table 2: Cumulative confirmed cases of SARS in Amoy Gardens, Block E, March 2003 (from Ching, et al. (0)). 
Predicted cases based on parameter estimation Methods 1 and 2. 



Parameter 


Dimension 


Method 1 


Fixed or Estimated 


Method 2 


Bound or Constraint 


So 


people 


788 


Fixed 


788 


So =788 


h 


people 


0.4 


Estimated 


0.0 


h + P) = 4 


Po 


people 


3.6 


Fixed by Iq 


4.0 


h + Po = 4 


(3 


days _1 x people -1 


1.4850/N 


Fixed 


0.8728/N 


< /? < 0.01 


b 




0.9315 


Estimated 


0.9994 


0.8 < b< 1 


V\ 


days -1 


0.9997 


Estimated 


0.2500 


vi = 0.25 


V2 


days -1 


0.7294 


Fixed by 6, fi, and Rq 


0.1000 


V 2 = 0.1 


Rq 




1.5154 


Fixed 


3.4765 


1.5 < Rq < 4 



Table 3: Estimated parameters for the SIPR model and the SARS outbreak in the Amoy Gardens, Hong Kong, 
March 2003, using Methods 1 and 2. Fixed parameters were found in Ching et al. ([2j). 



formula to estimate the percent of the population that 
is susceptible during an epidemic. The number who are 
susceptible can be computed using the Rq value and the 
percent of the population who became infected. Riley, et 
al. ([13]) determine the Rq value specifically for the SARS 
outbreak in Hong Kong is in the region 2.2 - 3.7. Using 
this, we compute that the number of susceptible individ- 
uals falls in the range 1.8 million - 3.1 million, where the 
smaller population correspond to larger R values. We 
assume a large Pn value, and, therefore we use the lower 
value, N = 1.8 million. 

For parameter estimation Method 1, we use the param- 
eters determined by fitting the data sets to the classical 
SIR model. The fixed and estimated parameters of the 
SIPR model are given in Tables [5] and [51 corresponding 
to time period one and three, respectively. The result- 
ing number of confirmed cases of SARS are presented in 
Table H 

For parameter estimation Method 2, we survey litera- 
ture to determine appropriate bounds on all of the SIPR 
model equations. We assume that the total number of 
infected individuals at t = is 1 for time period one, 
corresponding to the index case, and 470 for time period 
threqj, corresponding to the number of confirmed cases 
on March 29, 2003. The bounds and constraints for each 
parameter as well as the estimated parameter values are 
given in Tables [5] and [6] The resulting number of con- 
firmed cases of SARS are presented in Table 0J 



5 Some of these 470 confirmed cases would have already been re- 
moved from the disease transmission dynamics via recovery, death, 
or removal to isolation. 



3.3 Summary Results 

We summarize the infections for each of the data sets and 
their fitted parameters in Table [7] Specifically, the table 
contains the total number of individuals who became in- 
fected, the total number who became infected through 
contact with a regular infected (I) and super-spreading 
individuals (P), the total number of regular infected (J) 
and super-spreading individuals (P), and the individual 
reproduction numbers for the I and P classes (as SIR 
models). 



4 Discussion 

We have presented a modification of the classical SIR 
disease model that captures the effect of super-spreading 
individuals on an infectious disease epidemic. Using an 
idea from the progression of the SARS outbreak, we dis- 
tinguish the regular infected individuals from the super- 
spreading individuals by how long they remain outside 
of isolation; the super-spreading individuals spend longer 
outside of isolation than most infected individuals. The 
model was fit to data from the SARS epidemic, using two 
different parameter estimation methods. 

Parameter estimation Method 1, used parameters esti- 
mated for the classical SIR model. The resulting pa- 
rameters show only a slight super-spreading behavior in 
all cases studied. For Amoy Gardens and the short-term 
outbreak in Hong Kong, the / and P classes are not dis- 
tinguished by how long each stays outside of isolation. 
They are also not distinguished by their individual re- 
production numbers. In each case, the two recovery rates 
are similar to eachother, and similar to the recovery rate 



Day 


Feb 21 


Mar 17 


19 


21 


29 


Apr 12 


26 


May 10 


24 


Jun 6 


20 


Confirmed 


1 


95 


150 


203 


470 


1108 


1527 


1674 


1724 


1750 


1755 


SIR 


1 


100.9 


143.8 


204.9 


470.0 


1081.5 


1355.8 


1478.9 


1534.0 


1557.6 


1569.3 


Method 1 


1 


100.3 


142.9 


203.6 


470.0 


776.1 


983.6 


1124.8 


1221.0 


1282.6 


1328.4 


Method 2 


1 


103.8 


144.9 


202.2 


470.0 


1156.5 


1486.1 


1644.4 


1720.5 


1755.1 


1773.6 



Table 4: Cumulative confirmed cases of SARS in Hong Kong. Predicted cases based on parameter estimation for the 
classical SIR model, and the SI PR model using Methods 1 and 2. 



Parameter 


Dimension 


Method 1 


Fixed or Estimated 


Method 2 


Bound or Constraint 


So 


people 


1,799,999 


Fixed 


1,799,999 


5 =N-1 


h 


people 


1.0 


Estimated 


1.0 


h + Po = 1 


Po 


people 


0.0 


Fixed by Iq 


0.0 


Io + Po = 1 


P 


days -1 x people -1 


0.2586/N 


Fixed 


0.3738/N 


< /3 < 1/JV 


b 




0.6489 


Estimated 


0.8000 


0.8<6< 1 


V\ 


days -1 


0.0836 


Estimated 


0.2500 


v x = .25 


V2 


days -1 


0.0794 


Fixed by b, v\, and Ro 


0.1000 


z/ 2 = .1 


Ro 




3.1511 


Fixed 


1.9438 


1.5 < Ro < 4 



Table 5: Estimated parameters for the SI PR model and the SARS outbreak for February 21 
Hong Kong, using Methods 1 and 2. 



March 21, 2003, in 



Parameter 


Dimension 


Method 1 


Fixed or Estimated 


Method 2 


Bound or Constraint 


So 


people 


1,799,530 


Fixed 


1,799,530 


S =N-470 


Io 


people 


277.2 


Estimated 


311.6 


Io + Po = 470 


Po 


people 


192.8 


Fixed by Io 


158.4 


Io + Po = 470 


P 


days -1 x people -1 


0.1351/N 


Fixed 


0.1473/N 


< p < 1/JV 


b 




0.5237 


Estimated 


0.8910 


0.8 < 6 < 1 


V\ 


days -1 


25.6988 


Estimated 


0.2500 


v x = .25 


V2 


days -1 


0.0920 


Fixed by b, v\, and Ro 


0.1000 


Vl = .1 


Ro 




0.7025 


Fixed 


0.6851 


< Ro < 1 



Table 6: Estimated parameters for the SI PR model and the SARS outbreak for March 29 
Kong, using Methods 1 and 2. 



June 20, 2003 in Hong 





Amoy 


Gardens 




HK 


Period 1 




HK 


Period 3 






SIR 


M 1 


M 2 


SIR 


M 1 


M 2 


SIR 


M 1 


M 2 


Total Infected 


477.3 


479.1 


765.8 


1,709,770 


1,709,769 


1,405,410 


1578.9 


1426.1 


1791.5 


Contact with / 


473.3 


426.1 


742.1 


1,709,769 


1,097,524 


912,542 


1108.9 


4.1 


876.6 


Contact with P 




49.0 


19.7 




612,244 


492,867 




952.0 


444.9 


Total / 


477.3 


443.0 


761.4 


1,709,770 


1,109,516 


1,124,365 


1578.9 


777.9 


1489.2 


Total P 




36.1 


4.4 




600,253 


281,045 




648.2 


302.3 


Ro for / 


1.5154 


1.4779 


3.4735 


3.1511 


3.0944 


1.4953 


0.7025 


0.0053 


0.5889 


Ro for P 




2.0255 


8.6838 




3.2560 


3.7382 




1.4691 


1.4721 



Table 7: Summary of the outbreaks for each of the data sets and their fitted parameters: the total number of 
individuals who became infected; the total number who became infected through contact with a regular infected (/) 
and super-spreading individuals (P); the total number of regular infected (/) and super-spreading individuals (P); 
and the individual reproduction numbers for the / and P classes (as SIR models). 



of the SIR model estimated parameters. It is not sur- 
prising, therefore, that the final number of infected indi- 
viduals, given by the final number of individuals in the 
removed class, in the two models are within just a few 
individuals. 

The fitted parameters for the long-term outbreak in Hong 
Kong using parameter estimation Method 1 indicate that 
regular infected individuals move so quickly into isola- 
tion that we can disregard their influence on the disease 
spread. As evidence, over the entire course of the disease, 
they collectively only infect a total of 4.4 individuals. Due 
to the incredibly short time spent in the infected class, 
the individual reproductive number in this case is very 
low. 

Though the outbreak in the entire Hong Kong Special 
Administrative Region was split into three time periods 
corresponding to three different disease dynamics, the fit 
parameters do not match the third time period data set 
as well as one might hope. This is noticed for the classical 
SIR model, and, since the SIR parameters are used for 
estimation Method 1, we also see this in the fit parameters 
for Method 1. One possible fix for this problem is to 
split the third time period into other periods in which 
the outbreak has common dynamics. 

Parameter estimation Method 2, used research to set a 
priori bounds on the parameters of the SIPR model. 
In every case, the resulting parameters show super- 
spreading behavior, for example, the two infected classes 
are distinguished by their individual reproduction rates. 
All three cases show that, on average, each regular in- 
fected individual infects less than 1 other individual, 
while each super-spreading individual infects more than 
1 other individual. 

Considering the specific case of the outbreak at the Amoy 
Gardens apartment complex, Method 2 predicts that, 
without any other intervention, almost all of the residents 
of Block E, will become infected. In this case, both the 
regular infected and super-spreading individuals have ba- 
sic reproduction numbers larger than one. It is clear from 
the dire outcome predicted for the residents of the Amoy 
Gardens, Block E, super-spreading individuals must be 
brought into isolation as quickly as possible. 

For the long-term outbreak in Hong Kong (in the third 
time period) , the overall individual reproductive number 
is less than 1, which matches the notion that individu- 
als were taking precautions to protect themselves, and 
the disease spread was slowing down. Using the SIPR 
model, we see that the regular infected individuals do 
have an individual reproductive number less than 1, how- 
ever, the super-spreading individuals have a reproductive 
number larger than 1. Again, we see that it is impera- 
tive that the super-spreading individuals be brought into 
isolation as quickly as possible. 



Finally, it is evident from the short-term outbreak in 
Hong Kong (in the first time period), that without any 
control measures the spread of SARS in Hong Kong would 
have been extreme. The model predicts that more than 
1.7 million residents would have been infected over the 
course of the disease outbreak. Correspondingly, they 
would have experienced a large disease related mortality. 
The death rate for SARS is estimated to be around 10%, 
and so, Hong Kong would have lost an estimated 170,000 
residents. Noting that the final death toll in Hong Kong 
was (though tragic) only 299 people, the control measures 
put in place by the government of Hong Kong saved thou- 
sands of lives. 

The SIPR model is versatile; it can be used to examine 
an outbreak of any disease known to have super-spreading 
individuals, measles for example (see ( [T3])) . On the 
other hand, the model was built using the idea that 
super-spreading individuals stay out of isolation longer 
than regular infected individuals, as during the SARS 
epidemic. (There are documented cases of SARS infected 
individuals violating strict isolation mandate^.) For dis- 
eases where the super-spreading behavior is a result of 
differing transmission rates, one should use the Kemper 
model ([6]). 

Super-spreading individuals for infectious diseases pose 
a serious threat to public health. The SIPR model 
clearly demonstrates that infectious individuals must be 
removed from interactions with the susceptible individu- 
als as quickly as possible to decrease the seriousness of 
an infectious disease epidemic. 
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